The Meissner effect and massive particles as witnesses of macroscopic entanglement 
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We show how the spatial macroscopic entanglement equiv- 
alent to the off diagonal long range order (ODLRO) implies 
the Meissner effect and flux quantisation for a superconduc- 
tor. It is argued by analogy with superconductors that the 
Higgs field must also be entangled in the same way. Inter- 
nal (spin) entanglement is shown to be irrelevant within this 
context, although it can of course also be computed. 



Introduction. We have recently argued that high tem- 
perature macroscopic entanglement is possible and linked 
it to high T c superconductivity [1]. We have also dis- 
cussed the relationship between the notion of long range 
off diagonal order (ODLRO) [2] in a state and the exis- 
tence of bipartite and multipartite entanglement in the 
same state. Now we intend to extend this line of thought 
and show that the said multipartite entanglement implies 
two typical superconducting effects: the exclusion of the 
magnetic field from a simply connected superconductor 
(the Meissner effect) and the quantization of flux in mul- 
tiply connected regions. Although we will use the Hub- 
bard related models to aid the discussion in this paper, 
our results hold for any model which displays ODLRO, 
i.e. normal superconductors and superfluids. The paper 
will be finished with a speculation, maintaining that if 
Higgs bosons are found to exist, and if they are respon- 
sible for mass generation through the symmetry break- 
ing mechanism, then they must also be entangled. The 
reason is that the condensation of Higgs bosons is un- 
derstood to be the most likely mechanism for mass gen- 
eration in local gauge field theories (i.e. the Standard 
Model). 

Setting the scene. The model we analyse consists of a 
number of lattice sites, each of which can be occupied by 
fermions of spin up or spin down. Since fermions obey 
the Pauli exclusion principle, we can have at most two 
fermions attached to one and the same site. Let us in- 
troduce fermion creation and annihilation operators, cf s 
and a :S respectively, where the subscript i refers to the 
i-th lattice site and s refers for the value of the spin, | 
or j. The c operators satisfy the anticommutation rela- 



tions: {ci 7S ,Cj t } — 5ij5 s j, and c's and c^'s anticommute 
as usual. (Some general features of fermionic entangle- 
ment were analysed in [3-5]). 

We only need assume that our model has the interac- 
tion which favours formation of Cooper pairs of fermions 



of opposite spin at each site - these states are known as r\ 
states [6] and will be discussed below. The actual Hamil- 
tonian is not relevant for our present purposes. Suffice it 
to say that the r\ states are eigenstates of the Hubbard 
and realted models relevant for superconductivity [6,7]. 
Introducing r\ states. Suppose that there are n sites 
and suppose, further, that we introduce an operator 
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(1) 



that creates a coherent superposition of a Cooper pair in 
each of the lattice sites. This operator can be applied 
to the vacuum a number of times, each time creating 
a new coherent superposition. However, the number of 
applications, fc, cannot exceed the number of sites, n, 
since we cannot have more than one pair per site due to 
the exclusion principle. We now introduce the following 
basis 



|fc, n — k) := 



-1/2 



fa + ) fc |o> , 



(2) 



where the factor in front is just the necessary normal- 
isation. Here, the vacuum state |0) is annihilated by 
all c operators, Ci yS \0) = 0. We note that the originally 
defined rj operators can also have phase factors depen- 
dent on the location of the site on the lattice, like so 
Vk = J2 n elkn cli j c n j_- All the states generated with any 
r)k from the vacuum have the same amount of entangle- 
ment so that the extra phases will be ignored in the rest 
of the paper. However, the phase must be chosen, for 
otherwise, if we average over all possible phases, the re- 
sulting state is no longer entangled. Choosing a phase 
amounts to "symmetry breaking" and we will have more 
to say about it below. 

We can think of the rj states in the following way 
[1]. Suppose that k — 2. Then this means that we 
will be creating two repairs in total, but they cannot 
be created in the same lattice site. The state |2,n — 2) 
is therefore a symmetric superposition of all combina- 
tions of creating two pairs at two different sites. Let us, 
for the moment, use the label when the site is unoc- 
cupied and 1 when it is occupied. Then |2,n — 2) = 

(|00...11) + ...|11...000})/y^), i.e. the state is an equal 
superposition of states containing 2 states |1) and n — 2 
states |0). These states, due to their high degree of 
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symmetry, are much easier to handle than general arbi- 
trary superpositions and we can compute entanglement 
for them between any number of sites [8] . In this descrip- 
tion each site effectively holds one quantum bit, whose 
signifies that the site is empty and 1 signifies that the 
site is full. 

ODLRO. The main characteristic of rj states is the ex- 
istence of ODLRO, which implies its various supercon- 
ducting features, such as the Meissner effect and the flux 
quantisation [9]. The ODLRO is defined by the off di- 
agonal matrix elements of the two-site reduced density 
matrix being finite in the limit when the distance be- 
tween the sites diverges. Namely, 

lim (ct A ,c i;i c i;T ) — ► a (3) 

where a is a constant (independent of n) [2]. We will 
show that although the existence of off diagonal matrix 
elements does not guarantee the existence of entangle- 
ment between the two sites, it does guarantee the ex- 
istence of multi-site entanglement between all the sites. 
Note that here, by "correlations" we mean correlations 
between the number of electrons positioned at different 
sites i and j. This is different from spin-spin correla- 
tions, which would look at the occurrences of both elec- 
tron spins being up or down, or one being up and the 
other being down [4]. 

The r\ states are always of the form \k,n — k) := 

(S'|000...11))/ \J (2) , where S is the total symmctrisation 
operator. The reason why rj states are important for su- 
perconductivity is that this phenomenon can be under- 
stood to arise through a condensation of Cooper pairs. 
Condensation means that the temperature is so low that 
all particles are spread across the whole system - i.e. 
their wavelengths are as large as the system - so that 
all their wave functions overlap to a high degree (using 
the language of the first quantisation). This is why the 
r\ states are a good description as they represent equal 
superpositions of Cooper pairs across all the sites. 

We would now like to start to compute the entangle- 
ment between every two sites in the state \k,n — k). A 
simpler and more insightful task would be first to tell if 
and when every two qubits in a totally symmetric states 
are entangled. For this, we need only compute the re- 
duced two-qubit density matrix which can be written as: 

Pl2 {k) = a|00)(00| + 6|11)(11| + (#+)(V>+| (4) 

where - (|00) + |ll»/>/2 and a = 

^ (n-fc)(n-fc-i) _ 2k(n-k) ^ can easily check that 

n{n— 1) n[n— 1J J 

a + b + c= 1 and so the state is normalised. This density 
matrix is the same no matter how far the two sites are 
from each other, since the state is symmetric, and must 
therefore be identical for all qubits. We can easily test 
the Peres-Horodecki (partial transposition) condition for 



separability of this state [1]. This leads to states pi2(k) 
being entangled if and only if a+b— \J {a — b) 2 + Ac 2 < 0, 
which leads to (k — l)(n — k — 1) < k(n — k). This is, 
of course, satisfied for all n > 2 and 1 < k < n — 1. 
So, apart from the case when the total state is of the 
form |000..0) or |lll..l), there is always two-qubit en- 
tanglement present in symmetric states. Note, however, 
that in the limit of n and k becoming large no matter 
what their ratio may be - the value of the left hand side 
approaches the value of the right hand side and entan- 
glement thus disappears in the thermodynamical limit 

[!]■ 

Macroscopic (spatial) entanglement. The two 

point correlation function used in the calculation of the 
ODLRO in eq. (3) is, in fact, just one of the 15 numbers 
we need for the full two-site density matrix. In our sim- 
plified case of symmetric states in the ^-pairing model, 
this off diagonal clement is equal to c. However, for the 
density matrix we still need to know a and 6, and these 
numbers clearly affect the amount of entanglement. So, 
the first lesson is that two-site entanglement is not the 
same as the existence of ODLRO, and therefore two-site 
entanglement is not relevant for superconductivity. This 
does not mean, of course, that there is no entanglement 
in the whole of the lattice. In fact, the ODLRO implies 
that the two site density matrix contains classical cor- 
relations. This, together with the fact that the overall 
state of all electrons is pure, means that there is always 
bipartite entanglement present. For example, entangle- 
ment exists in the thermodynamical limit between two 
bunches of sites containing k and n — k sites respectively, 
and each bunch containing a site from our two site den- 
sity matrix [1]. 

The Meissner effect. Suppose that we exchange two 
electron pairs, one at the site 1 and one at the site n. 
We can imagine doing this adiabatically, although the 
requirement of adiabaticity is by no means necessary (it 
is merely convenient, as this evolution then generates no 
other effect apart from the one we wish to concentrate 
on). Suppose that whatever the total state is, the re- 
duced density matrix of sites 1 and n has a non- vanishing 
component of the state 

|*) = |0i)|l n ) + |li)|0„). (5) 

which, from the above discussion, means that we assume 
ODLRO. Then, after the swap, this component will look 
like 

|*) = |0 1 )|l„) + e i *|l 1 )|0„), (6) 

where $ = J Adl is the line integral of the vector poten- 
tial along the path traversed by the electron pair (with 
proper units introduced below). The reason why the two 
states in the superposition acquire different phases is that 
the electron pairs in two states undergo evolutions in op- 
posite directions of each other - this is, in fact, the well 
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known Aharonov-Bohm phase [10] . So, if in the first state 
the pair takes one path (i.e. the electron pair from site n 
moves to site 1), in the second state it takes the reverse 
of the same path (i.e. the electron pair from site 1 moves 
to site n). The off-diagonal element in the |0), |1) basis of 
the two site density matrix of sites 1 and n undergoes the 
following transformation c — ► e J *c. However, the overall 
state must be totally symmetric, and therefore 

e l<s> = 1. (7) 

From this we can conclude that 

* = f c l Adl = f c l l BdS = 2 ™- ® 

But, in the two (three) dimensional space, the electron 
pair can take any trajectory, and the only possible choice 
that satisfies the above is B = and n = 0. Therefore, 
in a connected region of a material exhibiting bipartite 
entanglement there is no magnetic field present. This is 
the Meissner effect. We note that it is well known that 
the magnetic field does penetrate the superconductor to 
a very small degree (falling off exponentially with the 
distance from the surface), but this cannot be explained 
with the present simple formalism and we need a more 
elaborate electrodynamic treatment. 
Flux Quantisation. Imagine instead that the region 
is not simply connected and that there is a hole in the 
middle pierced through by a magnetic field (there could 
be more than one hole and the same conclusion will hold 
for each of them) . Then its flux must be quantised. This 
follows immediately from eq. (8), since now not all paths 
are allowed. Namely, the field is now confined to the 
region where electrons cannot go, so that 

Therefore, we must have that ff 3> c = 2n7r, and so the 
flux is quantised in units of hc/2e: 

1>c = ng. (10) 

This is the flux quantisation effect. Note that the de- 
nominator contains twice the electron charge and this is 
a consequence of electrons forming Cooper pairs. 

The flux quantisation that we have just derived lies 
behind the persistent flow of electrical current in a su- 
perconductor. The flux is a consequence of the flowing 
current and any (continuous) dissipation cannot change 
the current continuously as the flux is discrete. Therefore 
the current persists indefinitely. 

Finally, if there is no ODLRO, meaning that as n — > oo 
we have that c — > in eq. (4) , then neither of the above 
two effects follow. Any phase now gained upon exchange 
of electrons as described above will not be reflected in 



the two site state and therefore we cannot argue that this 
phase has to have a special value. Therefore, bipartite 
entanglement is necessary for the Meissner effect and flux 
quantisation. Note briefly that the converse is not true. 
Not all entanglement will lead to superconductivity. For 
example, look at the state of two sites of the form: 1 00) + 
1 11). When we exchange the pairs we get no extra phase 
in the second ket (because they have opposite signs and so 
their product equals identity), so that the state remains 
the same. Therefore, any magnetic field is allowed to 
permeate such a state. This is why the ODLRO concerns 
the coherences between states 1 01) and |10). 
Mass from entanglement between Higgs bosons. 
We would now like to talk about the Higgs mechanism as 
the main explanation for the appearance of mass in local 
gauge field theories (see Weinberg [11] for a comprehen- 
sive introduction). In the modern field theory, gauge 
invariance of the Hamiltonian (or Lagrangian, which is 
more typically used) is invoked to explain the appear- 
ance of fields and their bosonic mediators. By "gauge 
transformation" we mean a transformation that acts on 
the wavefunction (or the field, more precisely) in the fol- 
lowing manner \^{x,t)) — ► e ^ x,t ^\^(x, t)), where 9{x,t) 
is just a phase that dependents both on space and time 
(i.e. it is local). In order for the Hamiltonian to remain 
invariant under this local change, we need to introduce 
an extra (vector) field, whose features exactly cancel out 
the effects of the local phase change. The necessary field 
turns out to be the electromegnetic field and its bosons 
are, of course, photons. The important point is that if we 
are to derive other forces from local gauge invariance (this 
requires phases that are non-commuting - i.e. matrices, 
but the concept is the same), the resulting bosons will 
always be massless. This result is intuitively clear: the 
local phase change has to be matched between arbitrary 
points and times and this therefore requires an infinite 
range force. The mediators of infinite range forces have 
to be massless. So, it appears that local gauge invariance 
cannot explain forces whose mediators are massive. 

A solution to this problem was found by Higgs [12] (and 
a number of other people, but Higgs was most promi- 
nent). The idea is that in addition to specifying the 
Hamiltonian of the fields, we also need to specify their 
actual physical state. This state need not possess the 
same symmetry properties of the Hamiltonian (hence the 
phrase "symmetry breaking" [11]). Suppose now that our 
local gauge invariance leads to several interacting mass- 
less fields. Suppose also, that one of the fields - known 
as the Higgs field - has condensed. In a mechanism that 
is completely analogous to the Meissner effect, the other 
fields will now be "expelled" from the Higgs field and 
will therefore become short range. In other words, their 
mediators will become massive. The condensation of the 
Higgs field therefore provides a mechanism to maintain 
local gauge invariance and have massive gauge fields at 
the same time, thereby circumventing previously men- 
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tioned limitation (the Higgs boson also acquires a mass 
in this process). Whether this is the correct way of ex- 
plaining the origin of mass in the Universe is still unclear 
as the search for Higgs bosons has so far been fruitless. 
However, one conclusion we can draw with more con- 
fidence (following this paper) is that if the Higgs field 
exists, then its bosons must be entangled. The reason 
is that the ODLRO, necessary for condensation, also im- 
plies existence of entanglement, and this is also true for 
Higgs condensation. The (obvious) fact that there are 
massive objects in the Universe would then be an entan- 
glement witness of the purely quantum correlations in 
the underlying Higgs field. 

Here we have to exercise some caution. The entangle- 
ment in the Higgs field is something that we refer to as 
"continuous variable entanglement" , as opposed to the 
discrete degrees of freedom of the Hubbard model above, 
and this quantity can become infinite. However, this is 
not a serious problem because with enough care this in- 
finity can always be controlled. 

We now show in a very simple example the connection 
between mass and entanglement that is meant to sub- 
stantiate the above dicussion, but is by no means a proof 
of it. Suppose we have a massive bosonic free field, (f> 7 
with the usual Lagrangian density l/2((d fl <fi) 2 + m 2 (f) 2 ) 
(this is in 1 + 1 dimensions), where m is the (fixed fi- 
nite) mass. This is an infinite continuous system and 
we divide it into two halfs (arbitrarily). Tracing one 
part out and computing the von Neuman entorpy of the 
remaining part results in the entropy of entanglement 
of E w lnl/m 2 a 2 [13], where a is some cutoff used to 
avoid the ultraviolet infinity (this divergence may also 
be avoided by using the relative entropy with respect 
to some coarse graining [14], much in the same way as 
Gibbs did in classical statistical mechanics) . The amount 
of entanglement clearly depends on the mass which could 
therefore be said to witness it. 

Spin Entanglement. We would like to finally point 
out an interesting curiosity that clarifies the notion of 
entanglement we have analysed in this paper. Namely, 
as we mentioned before, the relevant entanglement for su- 
perconductivity (and Higgs bosons) is the spatial entan- 
glement between numbers of electrons at different space 
points. What about the entanglement between the spin 
degrees of freedom? If the electrons occupy the same site, 
then they have to be anticorrelated (in the singlet state) 
because of Pauli's exclusion principle. If the electrons 
are on different sites, then they are not spin correlated 
in the r\ state (in the BCS model they would be, for a 
suffucicntly small distance [15]). This is because if we 
measure an electron in one site and then in another all 
four possibilities for their internal states are equally likely 
and so there can be no spin correlation present. There 
is a limit, however, in which the spin entanglement be- 
comes relevant. This is when the interaction between 
sites dominates the hopping amplitude in the Hubbard 



model and in the state where we have one electron per 
site. Here there is no ODLRO and the state is not super- 
conducting. However, the effective interaction between 
electrons at neighbouring sites is now of the Heisenberg 
type, as electrons can still exchange their locations and 
could at the same time have opposite spins. Therefore, 
at low temperatures there would be some two site spin 
entanglement present in the model, which is albeit not 
important for superconductivity. 

Conclusions. In one of our previous publications [1] 
we argued that macroscopic entanglement exists at high 
temperatures and is related to high temperature su- 
perconductivity. In the present work we showed that 
the consequences of that entanglement are the standard 
features of superconductors: the Meissner effect and 
flux quantisation. Therefore any experiments confirming 
these two effects are also automatically offering evidence 
for macroscopic entanglement. We have speculated that 
if the Higgs mechanism for mass generation is proven to 
be correct, then the resulting Higgs bosons will be found 
to be entangled. Be that as it may, one question remains 
open, both for superconductors, or for any other more 
general field. Can we extract this existing entanglement 
and use it for information processing? This would be 
very useful in practice, and it would seem that natural 
macroscopic entanglement could offer an infinite amount 
of quantum non-locality for genuine quantum informa- 
tion processing. This is the subject of an ongoing re- 
search. 
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